function [M,Mxp,Myp,Mxpx,Mxpxp,Mxpy,Mxpyp,Mypx,Mypxp,Mypy,Mypyp] = numDerivM_Real(params)
 
% Computing the steady 
[auxOut,errorMes] = DSGEmodel_ss(params);
params = auxOut.params;
xss = auxOut.xssTrans;
yss = auxOut.yssTrans;
 
%Unfold params
BETTA= params.BETTA;
B= params.B;
CHI= params.CHI;
CHI0= params.CHI0;
THETA= params.THETA;
DELTA= params.DELTA;
ALFA= params.ALFA;
PHI= params.PHI;
PHIzero= params.PHIzero;
ZI= params.ZI;
KAPAw= params.KAPAw;
ETA= params.ETA;
RHOR= params.RHOR;
PHIpai= params.PHIpai;
PHIpai_1= params.PHIpai_1;
PHIy= params.PHIy;
PHIy_1= params.PHIy_1;
PHIc= params.PHIc;
PHIc_1= params.PHIc_1;
PHIl= params.PHIl;
NU= params.NU;
U0= params.U0;
U0d= params.U0d;
RHOz= params.RHOz;
RHOd= params.RHOd;
RHOn= params.RHOn;
RHOp= params.RHOp;
RHOa= params.RHOa;
Kss= params.Kss;
OUTPUTss= params.OUTPUTss;
PAIss= params.PAIss;
MUZss= params.MUZss;
Css= params.Css;
AA= params.AA;
lss= params.lss;
Rss= params.Rss;
Wss= params.Wss;
 
% Current values in the model;
% x;
c_ba1 = xss(1);
muz_cu = xss(2);
d_cu = xss(3);
n_cu = xss(4);
paistar_cu = xss(5);
a_cu = xss(6);
% xp;
c_ba1p = xss(1);
muz_cup = xss(2);
d_cup = xss(3);
n_cup = xss(4);
paistar_cup = xss(5);
a_cup = xss(6);
 
% y
c_cu = yss(1);
r_cu = yss(2);
pai_cu = yss(3);
w_cu = yss(4);
l_cu = yss(5);
output_cu = yss(6);
mc_cu = yss(7);
evf_cu = yss(8);
vf_cu = yss(9);
p1_cu = yss(10);
p1Real_cu = yss(11);
dc_cu = yss(12);
% yp;
c_cup = yss(1);
r_cup = yss(2);
pai_cup = yss(3);
w_cup = yss(4);
l_cup = yss(5);
output_cup = yss(6);
mc_cup = yss(7);
evf_cup = yss(8);
vf_cup = yss(9);
p1_cup = yss(10);
p1Real_cup = yss(11);
dc_cup = yss(12);
 
% Setting the dimension for the matrices at first order
ny  = auxOut.ny;
nx  = auxOut.nx;
Mxp   = zeros(1,nx);
Myp   = zeros(1,ny);
 
%Setting dimension for the matrices at second order
Mxpx   = zeros(nx,nx);
Mxpxp  = zeros(nx,nx);
Mxpy   = zeros(nx,ny);
Mxpyp  = zeros(nx,ny);
 
Mypx   = zeros(ny,nx);
Mypxp  = zeros(ny,nx);
Mypy   = zeros(ny,ny);
Mypyp  = zeros(ny,ny);
 
% START DISPLAYING M
M(1,1)=(BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING M
 
% START DISPLAYING Mxp
Mxp(1,2)=- (BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxp(1,3)=(BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING Mxp
 
% START DISPLAYING Myp
Myp(1,1)=-(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Myp(1,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Myp
 
% START DISPLAYING Mxpx
Mxpx(2,1)=(B^2*BETTA*CHI^2*exp(-d_cu)*exp(-muz_cu)*exp(-muz_cup)*exp(c_ba1)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-vf_cup)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpx(3,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
Mxpx(2,2)=- (B^2*BETTA*CHI^2*exp(-d_cu)*exp(-muz_cu)*exp(-muz_cup)*exp(c_ba1)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-vf_cup)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpx(3,2)=(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
Mxpx(2,3)=(BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpx(3,3)=-(BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING Mxpx
 
% START DISPLAYING Mxpxp
Mxpxp(2,2)=(BETTA*exp(-d_cu)*exp(2*muz_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 2)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*(CHI0 - CHI*(CHI0 - 1) + 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (2*B*BETTA*CHI*exp(-d_cu)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (B^2*BETTA*CHI*exp(2*c_cu)*exp(-d_cu)*exp(-2*muz_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2) - (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA^2*ALFA*BETTA*exp(-d_cu)*exp(2*muz_cup)*exp(-2*vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(ALFA - 1)*(CHI - 1)^2*(CHI0 - 1)^2*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(2/(ALFA - 1))*exp(muz_cup)^(2*(CHI - 1)*(CHI0 - 1) + 2)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (2*AA*ALFA*BETTA*exp(-d_cu)*exp(2*muz_cup)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*BETTA*exp(-d_cu)*exp(2*muz_cup)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*((CHI - 1)*(CHI0 - 1) + 1)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 2)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (2*AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mxpxp(3,2)=- (BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpxp(2,3)=- (BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpxp(3,3)=(BETTA*exp(-d_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING Mxpxp
 
% START DISPLAYING Mxpy
Mxpy(2,1)=- (B*BETTA*CHI^2*exp(2*c_cu)*exp(-d_cu)*exp(-muz_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (BETTA*CHI*exp(-d_cu)*exp(c_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - (B*BETTA*CHI*exp(-d_cu)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (B^2*BETTA*CHI*exp(2*c_cu)*exp(-d_cu)*exp(-2*muz_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2) - (AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) - (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mxpy(3,1)=(BETTA*CHI*exp(-d_cu)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mxpy(2,8)=(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1)) + (AA^2*ALFA*BETTA*exp(-d_cu)*exp(-2*vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(1/(ALFA - 1))*exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpy(3,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mxpy
 
% START DISPLAYING Mxpyp
Mxpyp(2,1)=(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2) + (AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mxpyp(3,1)=-(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mxpyp(2,9)=(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1)) + (AA^2*ALFA*BETTA*exp(-d_cu)*exp(-2*vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(ALFA - 1)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(2/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mxpyp(3,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mxpyp
 
% START DISPLAYING Mypx
Mypx(1,1)=(B*BETTA*CHI^2*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mypx(9,1)=(AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-vf_cup)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mypx(1,2)=-(B*BETTA*CHI^2*exp(-d_cu)*exp(-muz_cu)*exp(c_ba1)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mypx(9,2)=-(AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-vf_cup)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mypx(1,3)=(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mypx(9,3)=(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mypx
 
% START DISPLAYING Mypxp
Mypxp(1,2)=(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2) + (AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mypxp(9,2)=(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(CHI0 - CHI*(CHI0 - 1))*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1)) + (AA^2*ALFA*BETTA*exp(-d_cu)*exp(-2*vf_cup)*exp(d_cup)*exp(muz_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(ALFA - 1)*(CHI - 1)*(CHI0 - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(2/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
Mypxp(1,3)=-(BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mypxp(9,3)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mypxp
 
% START DISPLAYING Mypy
Mypy(1,1)=- (BETTA*CHI^2*exp(-d_cu)*exp(c_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1) - (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(c_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2);
Mypy(9,1)=- (AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) - (AA*ALFA*B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-vf_cup)*exp(c_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mypy(1,8)=(AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mypy(9,8)=(AA^2*ALFA*BETTA*exp(-d_cu)*exp(-2*vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(1/(ALFA - 1))*exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^(2*(CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mypy
 
% START DISPLAYING Mypyp
Mypyp(1,1)=(BETTA*CHI*exp(2*c_cup)*exp(-d_cu)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(CHI + 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 2) - (BETTA*CHI*exp(-d_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
Mypyp(9,1)=(AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mypyp(1,9)=(AA*ALFA*BETTA*CHI*exp(-d_cu)*exp(-vf_cup)*exp(c_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
Mypyp(9,9)=(AA*ALFA*BETTA*exp(-d_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (AA^2*ALFA*BETTA*exp(-d_cu)*exp(-2*vf_cup)*exp(d_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(ALFA - 1)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 2))/(exp(evf_cu)^(2/(ALFA - 1))*exp(muz_cup)^(2*(CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
% END DISPLAYING Mypyp
 
end
